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method. Coherent and dissipative dynamics of the soUtons has been observed. 
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The realization of Bose-Einstein condensation (BEC) 
of weakly interacting atomic gases [Q stimulates strongly 
the exploration of nonlinear properties of matter waves. 
This supports the new field of nonlinear atom optics, e.g., 
four wave mixing in BEC's |j, as well as the study of 
various types of excitations. Of particular interest are 
macroscopically excited Bose condensed states, such as 
vortices and solitons. Vortices, well known from the stud- 
ies of liquid helium Q , have recently been observed in two 
component gaseous condensates [Q. 

Soliton-like solutions of the Gross-Pitaevskii equation 
are closely related to similar solutions in nonlinear optics 
describing the propagation of light pulses in optical fibres. 
Here, bright soliton solutions correspond to short pulses 
where the dispersion of the pulse is compensated by the 
self-phase modulation, i.e., the shape of the pulse does 
not change. Similarly, optical dark solitons correspond 
to intensity minima within a broad light pulse [H. 

In the case of nonlinear matter waves, bright solitons 
are only expected for an attractive interparticle interac- 
tion (s-wave scattering length a < 0) |^], whereas dark 
solitons, also called "kink-states", are expected to ex- 
ist for repulsive interactions (a > 0). Recent theoretical 
studies discuss the dynamics and stability of dark solitons 
[Q"0 ^^ ^^^^ ^^ concepts for their creation |lll|-|l3|l. 

Conceptually, solitons as particle-like objects provide 
a link of BEC physics to fiuid mechanics, nonlinear optics 
and fundamental particle physics. 

In this Letter we report on the experimental inves- 
tigation of dark solitons in cigar-shaped Bose-Einstein 
condensates in a dilute vapor of ^^Rb. Low lying ex- 
cited states are produced by imprinting a local phase 
onto the BEC wavefunction. By monitoring the evolu- 
tion of the density profile we study the successive dynam- 
ics of the wavefunction. The evolution of density minima 
travelling at a smaller velocity than the speed of sound 
in the trapped condensate is observed. By comparison 
to analytical and numerical solutions of the 3D Gross- 
Pitaevskii equation for our experimental conditions we 
identify these density minima to be moving dark solitons. 



In our experiment, a highly anisotropic confining po- 
tential leading to a strongly elongated shape of the con- 
densate allows us to be close to the (quasi) ID situation 
where dark solitons are expected to be dynamically sta- 
ble [||. Parallel to this work, soliton-like states in nearly 
spherical BEC's of ^^Na are investigated at NIST |l|. 

Dark solitons in matter waves are characterized by a 
local density minimum and a sharp phase gradient of 
the wavefunction at the position of the minimum (see 
Fig.|l|a,b). The shape of the soliton does not change. 
This is due to the balance between the repulsive interpar- 
ticle interaction trying to reduce the minimum and the 
phase gradient trying to enhance it. The macroscopic 
wavefunction of a dark soliton in a cylindrical harmonic 
trap forms a plane of minimum density (DS-plane) per- 
pendicular to the symmetry axis of the confining poten- 
tial. Thus, the corresponding density distribution shows 
a minimum at the DS-plane with a width of the order of 
the (local) correlation length. A dark soliton in a homo- 
geneous BEC of density no is described by the wavefunc- 
tion (see |l^ and references therein) 
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with the position Xk and velocity Vk of the DS-plane, the 
correlation length /q = (47rano)~^/^, and the speed of 
sound Cs 
ForT 



= y^ 4:71 anoh/m^ where m is the atom mass. 
in ID, dark solitons are stable. In this case, 
only solitons with zero velocity in the trap center do not 
move; otherwise they oscillate along the trap axis [gj. 
However, in 3D at finite T, dark solitons exhibit ther- 
modynamic and dynamical instabilities. The interaction 
of the soliton with the thermal cloud causes dissipation 
which accelerates the soliton. Ultimately, it reaches the 
speed of sound and disappears [|l^. The dynamical in- 
stability originates from the transfer of the (axial) soli- 
ton energy to the radial degrees of freedom and leads to 
the undulation of the DS-plane, and ultimately to the de- 



struction of the soliton. This instabihty is essentially sup- 
pressed for solitons in cigar-shaped traps with a strong 
radial confinement [||, such as in our experiment Q. 

As can be seen from Eq.(|l|), the local phase of the dark 
soliton wave function varies only in the vicinity of the DS- 
plane, x ~ x^, and is constant in the outer regions, with 
a phase difference A(/) between the parts left and right to 
the DS-plane (see, e.g., Fig.|l|b). 
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FIG. 1. Density distribution (a) and phase distribution (b) 
of a dark soliton state with A(/) = tt. The density minimum 
has a width ^ Iq. The scheme for the generation of dark 
solitons by phase imprinting is shown in (c), where le is the 
width of the potential edge. 

To generate dark solitons we apply the method of phase 
imprinting |13], which allows one also to create vortices 
and other textures in EEC's. We apply a homogeneous 
potential Uint^ generated by the dipole potential of a far 
detuned laser beam, to one half of the condensate wave- 
function (Fig.[l|c). The potential is pulsed on for a time 
tp, such that the wavefunction locally acquires an ad- 
ditional phase factor e"*^"^, with A^ = Uinttp/h ~ tt. 
The pulse duration is chosen to be short compared to 
the correlation time of the condensate, tc = ^//i, where 
jii is the chemical potential. This ensures that the effect 
of the light pulse is mainly a change of the phase of the 
BEC, whereas changes of the density during this time can 
be neglected. Note, however, that due to the imprinted 
phase, at larger times one expects an adjustment of the 
phase and density distribution in the condensate. This 
will lead to the formation of a dark soliton and also to 
additional structures as discussed below. 

In our experimental setup (see Q), condensates con- 
taining typically 1.5 x 10^ atoms in the (F=2, 7tif=+2)- 
state, with less than 10% of the atoms being in the ther- 
mal cloud, are produced every 20s. The fundamental fre- 
quencies of our static magnetic trap are a;^^ = 27r x 14Hz 
and u;± = 27r x 425Hz along the axial and radial di- 
rections, respectively. The condensates are cigar-shaped 
with the long axis (x-axis) oriented horizontally. 

For the phase imprinting potential Uint^ a blue de- 
tuned, far off resonant laser field (A =532nm) of inten- 
sity / ^ 20W/mm^ pulsed for a time tp = 20/is results 
in a phase shift A0 of the order of tt p!?!] . Spontaneous 



processes can be totally neglected. A high quality optical 
system is used to image an intensity profile onto the BEC, 
nearly corresponding to a step function with a width of 
the edge, /g, smaller than 3/im (see Fig.|l|c). The corre- 
sponding potential gradient leads to a force transferring 
momentum locally to the wave function and supporting 
the creation of a density minimum at the position of the 
DS-plane for the dark soliton. Note that also the velocity 
of the soliton depends on /g (see Fig.||c). 

After applying the dipole potential we let the atoms 
evolve within the magnetic trap for a variable time tev 
We then release the BEC from the trap (switched off 
within 200/is) and take an absorption image of the den- 
sity distribution after a time-of-flight troF = 4ms (re- 
ducing the density in order to get a good signal-to-noise 
ratio in the images). 

In series of measurements we have studied the creation 
and successive dynamics of dark solitons as a function of 
the evolution time and the imprinted phase. Fig.|| shows 
density profiles of the atomic clouds for different evolu- 
tion times in the magnetic trap, tev The potential Uint 
has been applied to the part of the BEC with x < 0. For 
this measurement the potential strength was estimated 
to correspond to a phase shift of ~ tt. 

For short evolution times the density profile of the BEC 
shows a pronounced minimum (contrast about 40%). Af- 
ter a time of typically tev ^ 1.5ms a second minimum 
appears. Both minima (contrast about 20% each) travel 
in opposite directions and in general with different ve- 
locities. Fig.Pa) shows the evolution of these two minima 
in comparison to theoretical results obtained numerically 
from the 3D Gross-Pitaevskii equation. 

One of the most important results of this work is that 
both structures move with velocities which are smaller 
than the speed of sound (cg ~ 3.7mm/s for our parame- 
ters) and depend on the applied phase shift. Therefore, 
the observed structures are different from sound waves in 
a condensate as first observed at MIT [Q. We identify 
the minimum moving slowly in the negative x-direction 
to be the DS-plane of a dark soliton. 




tey= 0.5ms 2ms 3ms 4ms 9ms 12ms 

FIG. 2. Absorption images of EEC's with kink- wise struc- 
tures propagating in the direction of the long condensate 
axis, for different evolution times in the magnetic trap, tev 

(A0 ^7T, N ^1.5x 10^ tTOF = 4ms). 
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FIG. 3. (a) Position of the experimentally observed density 
minima (see Fig.g|) versus evolution time in the magnetic trap 
for A(/) ^ 0.57r. The dashed lines show results from the 3D 
simulation for A^ = 5 x 10^, le = 3/im, and Acj) = 27r/3. (b) 
and (c) show the dependence of Vk/cs on the imprinted phase 
(for le = 3/im) and on le (for Acj) = tt) obtained numerically 
in quasi ID simulations (see [16]) (with A^ = 1.5 x 10^). 

We have performed series of measurements with dif- 
ferent parameter sets for /g and the product of laser in- 
tensity and imprinting time. The velocity of the dark 
soliton could thereby be varied between v^ = 2.0mm/s 
(Fig.||a) and Vk = 3.0mm/s. For fixed /e, the velocity Vk 
decreases with increasing Acj). An increase of A0 ~ 7r/2 
by a factor of 1.5 results in decreasing Vk by 10%, in 
agreement with theoretical results (see Fig.||b,c). For 
significantly reduced imprinted phase we did not observe 
any dark soliton structures. For higher imprinted phase 
values more complex structures with several density min- 
ima were observed. 

In addition to the dark soliton, the dipole potential cre- 
ates a density wave travelling in the positive x-direction 
with a velocity close to Cg. After opening the trap, a 
complex dynamics results in the appearance of a second 
minimum behind the density wave as explained below. 

To understand the generation of dark solitons and their 
behavior in the initial stages of the evolution we have per- 
formed numerical simulations of the 3D Gross-Pitaevskii 
equation. Computing time limitations have restricted our 
studies to atom numbers below 5 x 10^. The simulations 
describe well the case T = 0, but ignore the effects of 
thermodynamic instability. The latter was analyzed by 
using a generalization of the theory of Ref. [|lO[ . 

Our theoretical findings are summarized as follows: 
1. The results of the simulations agree well with the ex- 
perimental observations. After applying a phase chang- 
ing potential, a dark soliton moves in the negative x- 
direction (Fig.^). The generation of the soliton by the 
phase imprinting method is accompanied by a density 
wave moving in the opposite direction. The maximum 
of the density wave travels with a velocity ^ Cg, inde- 
pendently of the values of Acj) and le. A characteristic 
time for the creation of the soliton is of order Iq/cs and 
in our case it does not exceed fractions of ms. Note that 
after this time the soliton-related phase slip in the wave- 
function is affected by a complex dynamics of the soliton 
generation and will be different from the imprinted Acj). 



2. For a fixed /g? the increase of A^ from tt to 27r, 37r, . . . 
leads to the creation of double, triple solitons etc. EEC's 
with several dark solitons were also observed experimen- 
tally and are currently investigated in detail. 

3. The initial soliton velocity decreases as /g ^ (Fig. 
^). As observed experimentally, typical soliton veloci- 
ties (for le < 3/im) are smaller than Cs and grow with 
the number of atoms. Velocities of the accompanying 
density waves are close to Cg. These waves move away 
from the center of the trap, broaden and eventually van- 
ish (Fig.0). This is in contrast to solitons, which are 
expected to oscillate in the trap, retaining their width 
and absolute depth. However, the observation of these 
oscillations would require longer lifetimes of the solitons 
(limited by dissipation to ~ 10ms, see below). Within 
this time scale we find no signatures of dynamical insta- 
bility and only reveal a moderate change (< 10%) of the 
soliton velocity, in agreement with our experiments. 

4. The situation changes after opening the trap and 
allowing the condensate to ballistically expand in the 
radial direction. The simulation shows that the soli- 
ton velocity decreases rapidly, while its width grows. 
To understand this aspect we have used a scaling ap- 
proach, similar to that of [|l9|,|0[, for the radial ballis- 
tic expansion of an infinitely elongated condensate con- 
taining a moving kink. This approach (valid for cjj^ < 
troF ^ m/^^1) predicts a soliton velocity VkitTOp) ^ 
Vfe(^e^)ln(2cj^tTOF)/^±^TOF, whcrc Vk{tev) IS the soli- 
ton velocity at tev just before switching off the trap. This 
result agrees very well with both experimental data and 
numerical simulations for a finite size BEC. Moreover, in 
a short time It of — ^^ after switching off the trap the 
density develops a second minimum located between the 
density wave and the dark soliton. This minimum has a 
width and depth comparable to those of the dark soliton 
(see Fig.Qb). Its position regarded as a function of tev 
moves with a constant velocity similar to that of the soli- 
ton. The creation of this second minimum is a coherent 
phenomenon and can be attributed to a dynamically ac- 
quired phase of the wave function in the region between 
the density wave and the dark soliton. 

The results of the experiment also show clear sig- 
natures of the presence of dissipation originating from 
the interaction of the soliton with the thermal cloud. 
We observe a decrease of the contrast of the soliton by 
~50% on a time scale of 10ms. This is in contradic- 
tion with the nondissipative dynamics, where the con- 
trast should even increase for the soliton moving away 
from the trap center. The soliton energy is then propor- 
tional to Uq (x/c)(l-'i;^/c^(x/c))^/^ (see |lo|) and should 
remain constant. This leads to a constant absolute depth 
of the soliton and hence gives the contrast proportional to 
(1 — v'^/c^g{xk)) ^ nQ{xk)~^ ' The decrease of the soliton 
contrast can therefore only be explained by the presence 
of dissipation decreasing the soliton energy. As the life 



time of the soliton is sensitive to the gas temperature, 
the studies of dissipative dynamics of solitons will offer a 
unique possibility for thermometry of EEC's in the con- 
ditions where the thermal cloud is not discernible. 




20 40 

x[|im] 



20 40 

x[|im] 



FIG. 4. Evolution of the density distribution obtained 
numerically from the 3D simulations for N = 10^ and 
Ac/) = 27r/3. (a) Evolution inside the magnetic trap for dif- 
ferent tev] the dark soliton is marked by an arrow, (b) After 
troF = 4ms, a second density minimum is observed behind 
the density wave. 

In conclusion, we have created dark solitons by a phase 
imprinting method and studied their dynamics. A de- 
tailed comparison to theory and numerical simulations 
allows us to identify the creation of dark solitons trav- 
elling with approximately constant velocity smaller than 
the speed of sound. The initial stages of the evolution 
and the radial ballistic expansion of the sample are well 
described by the T = approach which also shows the 
absence of dynamical instabilities. The decrease of the 
soliton contrast gives a clear signature of dissipation in 
the soliton dynamics. 

For the study of dark solitons with even lower veloc- 
ities, an initial density preparation of the BEC in the 
magnetic trap may be useful. This can be done, e.g., by 
applying adiabatically an additional blue detuned laser 
beam prior to the phase imprinting pulse. A promising 
attempt will be the realization of dark solitons in elon- 
gated dipole traps, e.g., generated by a blue detuned hol- 
low laser beam [^. With spin as an additional degree 
of freedom, the dynamics of dark solitons in condensates 
containing spin-domaines or spin waves can be studied. 

We expect that the study of soliton structures in EEC's 
opens a new direction in atomic physics, related to non- 
linear phenomena in a dissipative environment. 
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